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Density-functional Theory (DFT) approaches have recently been used to judge the topological 
order of various materials despite its well-known band gap underestimation. Use of the more accurate 
quasi-particle GW approach reveals here few cases where DFT identifications are false-positive, 
possibly misguiding experimental searches of materials that are topological insulators (TI) in DFT 
but not expected to be TI in reality. We also present the case of false-positive due to the incorrect 
choice of crystal structures and adress the relevancy of such choice of crystal structure with respect 
to the ground state one and thermodynamical instability with respect to binary competing phases. 
We conclude that it is then necessary to consider both the correct ground state crystal structure 
and the correct Hamiltonian in order to predict new TI. 

PACS numbers: 73.43.-f, 72.25.Hg, 73.20.-r, 85.75.-d 



3D bulk topological insulators (TI) have a band inver- 
sion, e.g. the s-likc conduction band below the (p, (i)-likc 
valence band at the time reversal invariant momentum 
such as r. This order of bands is defined through the in- 
version energy = Eg — £p,d being negative (see Figure 
1). Specifically, for normal insulators such as CdTe, 
is positive while for topological insulators such as HgTe, 
Ai is negative. The excitation gap defined as the en- 
ergy difference between the highest occupied state and 
the lowest unoccupied state can be either zero (Fig lb) 
or positive (Fig Ic) even in TI. When one tuncates a 3D 
bulk TI to create a 2D surface or interface, new states 
will appear inside the 2D excitation band gap: these will 
have linear dispersion ('massless fermions') crossing each 
other ('gapless') within the 2D excitation gap and main- 
tain spin-polarization without net charge current. 

Interest in their exotic properties has motivated recent 
search and discovery efforts of materials that would be 
topological insulators Theory predicts interesting 

cases expected to be TI's and pleas are then made to ex- 
perimentalists to make and probe such materials in the 
laboratory. As exciting as such a close level of theory- 
experiment interaction is, there are two potential issues 
that can impede progress. First, most of these calcula- 
tions rely on density functional theory (DFT) or variants 
of it these methods are known to systemati- 

cally underestimate band gaps specifically by placing the 
conduction band too low and the valence band too high, 
thereby creating the potential of 'false positive' predic- 
tion of a material being TI in DFT but not in reality, 
to the detriment of experiment-theory interaction. Sec- 
ond, sometimes [iHll theoretical predictions of negative 
inversion energy A^ are conducted on assumed crystal 
structures of given candidate materials; in some cases 
those assumed structures are energetically quite far from 
the ground state structure. This results in predicting an 



'Electronic address: | alex.zunger@gmail.com"] 



exciting property (e.g., TI) in unrealizable structures: we 
present here a case where the assumed crystal structure 
makes the compound thermodynamically unstable with 
respect to dissociation in competing binary phases. 

In this Communication, we revisit the topological or- 
der of several families of materials previously predicted 
to be TI within DFT, but now we are using (i) im- 
proved approximations for the electronic levels GW and 
(ii) insist on calculating topological order in structures 
that arc ground states or very close to it. The group 
of materials we study are: (a) binary compounds InAs 
and HgTe (b) Half Heusler compounds ABX [A=(Y,Lu); 
B=(Pt,Pd) and X=(Sb,Bi)] (c) Antiperovskite ni- 

tride (M)3BiN [(M=Ca,Sr,Ba)] 0, and (d) Honeycomb- 
lattice chalcogenides LiAgSe and NaAgSe 0- We first 
describe the method we used based on quasiparticle self- 
consistent GW. Then, we present a list of false-positive 
cases due to DFT and as false-positives due to the wrong 
crystal structure. Finally, we assess the potential of the 
recently used @ screened Hybrid functional, HSE06 
with respect to both DFT and GW methods. 

Methods: Contrary to DFT, quasi-particle (QP) ener- 
gies can be compared directly with experiments. Also, 
GW does not suffer from self-interaction error and is or- 
bital dependent [ll| , which are essential for excitation gap 
prediction. Out of the many fiavors of GW available, we 
employ the fully self-consistent (sc-) version of GW|12|. 
The need for self-consistence is obvious in the case of 
topological insulator: usually GW is employed within the 
so-called GqWo approximation, which is a perturbative 
correction to energies obtained from a mean-field Hamil- 
tonian (usually DFT). However, such approximation as- 
sumes implicitly that the mean-field Hamiltonian could 
describe the system of interest accurately enough so that 
its energies are considered to be good approximation to 
quasiparticle energies [l^. DFT dramatic underestima- 
tion of band gaps could result in metallic band structure 
within DFT for narrow band gap semiconductors. There- 
fore, Go Wo is not effective in such cases and one need to 
self-consistently update both QP energies and wavefunc- 
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tions in the GW Hamiltonian. Furthermore, spin orbit 
splittings calculated using norm-conserving relativistic 
pseudopotentials Q were added to the GW quasiparti- 
cle energies as a first order perturbation. All GW calcula- 
tions were performed using the crystal structure relaxed 
within generalized gradient approximation (GGA). 

False-positive predictions of topological characteristics 
due to DFT: Calculated DFT and sc-GW inversion en- 
ergies for previously proposed TI materials are given in 
Table I, together with their topological designation. 

(a) InAs and HgTe: InAs is a narrow-excitation-gap 
normal semiconductor with Af^^*^'= 0.43 eV. Sc-GW re- 
produces it well (see Table I) while DFT predicts a neg- 
ative inversion energy and consequently identify it incor- 
rectly as TI. 

Therefore, the assignement of topological order to InAs 
within DFT is considered a false positive and is due to 
the strong underestimation of the band gap in DFT. In 
the case of HgTe, both Af and Af^^ are negative, 
making HgTe a topological insulator in agreement with 
experimental findings. Only sc-GW is able to qualita- 
tively and quantitatively reproduce experimental inver- 
sion energy and excitation gap for both HgTe and InAs. 

(b) The Filled Tetrahedral Structures (FTS): ABX 
with A=(Y,Lu); B=(Pd,Pt) and X=(Sb,Bi) crystallize 
in a structure based on the zinc blende lattice which has 
a fourfold coordinated cation at position Ti = (0,0,0) and 
a fourfold coordinated anion at T2 = {j,j,j)- Whereas 
in the zinc blende lattice, the two tetrahedral intersitial 
sites - ITi at (^,^,5) coordinated by four T2 anions, 
and IT2 at (|, |, |) coordinated by four Ti cations - arc 
empty, the FTS are based on stuffing an additional atom 
either at ITi (called a-type FTS or at IT2 (/3- 
typeFTS We study both types. We found twice 
more false positive in a-type than in /3-type. Such dif- 
ference can be attributed to the localization of the VBM 
wavefunction: in a-type, VBM is formed by non -bonding 
states, highly localized around atoms while in /3-typc, 
VBM is formed by more dclocalized bonding states. As 
a result, a-type FTS suffers more from DFT-induced self- 
interaction error (T^]: VBM is predicted to lie too high 
in energy within DFT, closing the band gap and in some 
extreme cases inverting the order of the bands around 
Fermi level and therefore incorrectly assign it to be TI. 

(c) Strained antiperovskite nitrides M^BiN (M=Ca,Sr 
and Ba): The materials belonging to this family are 
said to be topological insulators under strain (but not as 
strain free). According to DFT calculations, unstrained 
MsBiN (M=Ca, Sr and Ba) antiperovskites have a small 
energy gap between low-lying unoccupied s-states of Bi 
and occupied p-states of Bi (defined as the inversion en- 
ergy Ai), prompting the suggestion Q that application 
of external strain can drive such small gap systems to 
have Tl-like negative inversion energy. Band gaps ob- 
tained within sc-GW for unstrained bulk are in good 
agreement with experimental data when available |18| . 
Such GW calculation for MgBiN (M=Ca,Sr,Ba) showed 
that the energy separation between p-likc valence states 




FIG. 1: (color online). Schematic of the band ordering at the 
vicinity of Fermi level for (a) normal insulator, (b) topological 
semimetal and (c) topological insulator. Red color represents 
(p,d)-like character while blue color represents s-like charac- 
ter. 



and s-like conduction states significantly increases from 
AfFT = 0.65 eV to Af"^^ = 1.85 eV for CagBiN, from 
AfFT ^ 0.26 eV to Af^^ = 1.55 eV for SrgBiN and 
from AfPT ^ i_Q^ ^sc-gw ^ ^ BagBiN, 

making gap closure by strain extremely unlikely. There- 
fore, strained (Ca,Sr,Ba)3BiN are probably /aZse positive. 

False Positive due to the wrong choice of crystal struc- 
ture: LiAgSe and NaAgSe in the assumed honeycomb- 
lattice (ZrBeSi-type structure see Fig 2a) are one of 
the few DFT-predicted TI materials with light elements 
{Z < 50). The VBM is formed by p-d hybridized states 
from Se and Ag while the conduction band minimum 
(CBM) is exclusively s-like. Within DFT, LiAgSe in 
the ZrBeSi-structure has a negative inversion energy at 
F while NaAgSe in the same structure presents an odd 
number of negative inversion energy at F, A and L, which 
makes both of them TIs within such crystal structure 
and DFT. Within sc-GW, none of the negative inver- 
sion energy are observed, leading to the conclusion that 
NaAgSe and LiAgSe are normal insulators in the ZrBeSi- 
structure (sec table I for values of inversion energy at F). 
This is similar to the false positive due to DFT presented 
in the previous section. 

Unfortunately, ZrBeSi-type structure assumed for 
(Li,Na)AgSe in reference is not the ground state struc- 
ture. According to our total energy calculations per- 
formed on more than 50 potential crystal structures, 
LiAgSe crystallizes in a LiCaN-structure (Figure 2b) 
(0.085 eV/atom lower than ZrBeSi-type) and NaAgSe 
crystallizes in the PbCl2-type structure (Figure 2c) (0.05 
eV/atom lower than ZrBeSi-type). Moreover, LiAgSe in 
the ZrBeSi-type structure is unstable with respect to dis- 
sociation into Li2Se and Ag2Sc. 

Interestingly, both DFT and sc-GW predict a posi- 
tive inversion energy for (Li,Na)AgSe in their respective 
ground state structure at any k-point in the Brillouin 
zone. Therefore, we discover a very special case of false 
positive due to the incorrect choice of ground-state struc- 
ture, independently of the type of methods used to pre- 
dict the inversion energies. 
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(a) 



(b) 



(c) 



ZrBeSi-type 
P63/mmc 



LiCaN-type 
Pnma 



PbClj-type 
Pnma 






FIG. 2: (color online). Structure type for (a) ZrBeSi-type (b) 
LiCaN-type and (c) PbCl2-type. 



Is HSE06 the answer ? 

Recently, Chen et al. Q determine new topological in- 
sulators based on screened hybrid functional HSE06 [l^ , 
which is known to improve band gaps of solids with re- 
spect to DFT [H]. We performed HSE06 calculations 
for aU the materials in Table I and CuTl(S,Se)2 [H 
recently predicted to be TI P. For most compounds, 
HSE06 gives qualitatively and quantitatively very similar 
results compared to sc-GW, with a 20% deviation for A; 
with respect to sc-GW results. However, two problematic 
cases appear where HSE06 is in qualitative disagreement 



with both DFT and sc-GW: HgTe where Af SEoe ^ q 
and the determination of topological order is then uncer- 
tain and CuTlSs where AfSEoe > q (see Table II) yet 
l^sc-GW ^ 0. Such false negative cases occurs because 
of the dependence of HSE06 on a fixed mixing parameter 
for exchange, Qx = j- The mixing parameter could 
be linked to inverse dielectric constant in the long range 
limit flx « ^ HH-ill. In the case of HSE06 and HgTe, 
a, = i and'elfTc ^ 2O.O. 

In conclusion, we revisited some of the numerous the- 
oretical predictions of TI and found that many of them 
are subject to false positive and false negative cases due 
to an inadequate use of ab initio methods or assumed 
structure. We use a robust method based on GW in or- 
der to avoid false positive and false negative. We adress 
the importance of using both the ground state structure 
and the correct Hamiltonian in order to search for new 
TIs. The failure to do so could result in false prediction 
and might mislead experimentalists in their effort to grow 
new TIs. 
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NREL. 
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TABLE I: Calculated DFT and GW inversion energies Ai, and topological class (T=topological and N=normal) for InAs, 
HgTe, (a-,/3-)(Y,Lu)(Pd,Pt)(Sb,Bi), (Li,Na)AgSe in the assumed ZrBeSi-type structure and their ground state structure. 



Inversion energy 

Topological ? 

Material (Structure) [gVl 





DFT 


olj V. I V V 


DFT 


DC Vjr VV 




-0.6 


0.65 
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N 


rlgie I^ZDj 


-i.i 




1 


1 




-0.21 


1.04 


T 


N 


YPdRi ('/^-FTS'l" 


-0.09 


0.5 


T 


N 


T,iiPtSh ('iT-FTS') 


-0 fi^ 


45 


T 


N 


LuPdBi (a-FTS) 


-0.27 


0.52 


T 


N 


YPtSb (a-FTS) 


-0.5 


0.45 


T 


N 


YPdBi (a-FTS) 


-0.1 


0.45 


T 


N 


LiAgSe (Fig. 2a)'' 


-0.17 


1.3 


T 


N 


LiAgSc (Fig. 2b) 


0.4 


1.9 


N 


N 


NaAgSe (Fig. 2a) 


-0.23 


0.8 


T 


N 


NaAgSe (Fig. 2c) " 


-0.0 


0.75 


N 


N 



"calculated ground state 
"calculated ground state 

''Excited state structure 0.085 eV/atom higher 
'^Ground state structure 

■^Excited state structure 0.05 eV/atom higher 
'^Ground state structure 



TABLE II: Calculated DFT, HSE06 and GW inversion energies Ai, and topological class (T=topological and N=normal) for 
CuTl(S,Se)2. CH stands for chalcopyrite. 





Inversion energy [eV] 
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DFT HSE06 sc 


-GW 


CuTlSz (CH) 


-0.77 0.06 -0.41 


T N 


T 


CuTlSea (CH) 


-1.0 -0.4 -0.85 


T T 
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